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Q_|' 1 INTRODUCTION 

^ ■ A bursting radio source, GCRT J1745-3009, was discovered 
' at 0.33 GHz in a radio monitoring program of the Galactic 

(^3 center region made on September 30, 2002 (Hyman et al. 

r* ' 2005a). Five ~ 10 min bursts with peak flux of ~ 1.67 Jy 
. were detected at an apparently regular period of ~ 77 min 
o from the source. Activity (only one single ~ 0.5 Jy burst) 

J_3 had been detected again by GMRT in 330 MHz at 2003 

Ctf \ September 28 (Hyman et al. 2005b). The source appears to 
be transient because it was not active at the 1998 Septem- 
ber 25 and 26 epochs of VLA observation, and had not been 
detected in some other epochs of observation in 2002 and 
2003. Observations indicate that (Hyman et al. 2005b) the 
burst detected in 2003 is an isolate one although additional 
undetected bursts occurred with 77 min period like the 2002 
bursts can not be completely ruled out. Assuming that the 
2003 burst is an isolated one, Hyman et al. (2005b) esti- 
mated crudely that the duty cycle of the transient behavior 
is about 10%. 

Given that (1) the source's brightness temperature 
would exceed 10 12 kelvin if it is farther than 100 pc; (2) the 
source's observational properties are not directly compatible 
with that of any known coherent emitters like white dwarfs 
or pulsars; Hyman et al. (2005a) concluded that it it is not 
likely to be a incoherent emitter but rather might be one of 
a new class of coherent emitter. Kulkarni & Phinney (2005) 
argued that the source could be a nulling radio pulsar, like 
PSR J 1752+2359 which has quasi periodic nulling behav- 
ior (Lewandowski et al. 2004). It is pointed out (Turolla, 



ABSTRACT 

A unique transient bursting radio source, GCRT J1745-3009, has been discovered (Hy- 
man et al. 2005a) near the direction of the Galactic center. It is still an open question 
to explain this phenomenon, although some efforts to understand its nature have been 
made. This paper shows that most of the observed features can be reproduced by our 
proposed precessing pulsar model. It is found that the precession angle of the pulsar 
should be larger (> 15°) than that of previously known precessing pulsars, which have 
a precession angle < 10°, if the beam width of the pulsar is larger than 10 degree. 
The pulsar could be a nulling (or even extremely) radio pulsars to account for the 
transient nature of the source. This model can be confirmed if a pulsar is detected 
at the position of the source. The pulsar could hardly be a normal neutron star (but 
probably a solid quark star) if the spin period of the pulsar is detected to be > 10 ms 
in the future. 
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Possenti & Treves 2005) that the phenomenon is compati- 
ble with what is expected from the interaction of wind and 
magnetosphere of two pulsars in a binary system. This sce- 
nario predicts: (1) a pulsar should be detectable at frequency 
higher than 1 GHz; (2) the X-ray luminosity from the shock 
should be 10 32 ergs s -1 which is too low to be detectable by 
contemporary facilities. The source could be a white dwarf 
(Zhang & Gil 2005), which may actually behave like a pul- 
sar and create the activity observed. This scenario predicts 
that deep IR exposure with large telescope may lead to the 
discovery of the counterpart of GCRT J 1745-3009. A conclu- 
sive understanding, however, has not been achieved yet, and 
could only be accomplished through further observation. 

An alternative effort is tried in this paper to explain the 
observational features of GCRT 1745-3009. We propose that 
the source could simply be a spinning pulsar precessing with 
a period of ~ 77 min. The duration and period of the bursts 
can be explained with a broad choice of parameters, as long 
as the precession angle is not very small (> 15 degrees). It 
is worth noting that the wobble angle of the pulsar could 
be typically of tens of degrees (Melatos 2000) if the free 
precession period is close to the radiation-driven precession 
period. Given that the brightness temperature could be as 
high as 10 28 — 10 30 K, a pulsar could reproduce the observed 
flux even if it is as far as 10 kpc away. The transient nature 
of the source would be understandable if the pulsar is an 
extremely nulling radio pulsar (Backer 1970; Ritchings 1976; 
Manchester). Some of the discovered nulling pulsars could 
have a huge nulling fraction. PSR 0826-34 is a case in point, 
whose nulling fraction is 70 ± 35 percent (Biggs 1992). PSR 
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B1931+24 switches off for ~ 90% of time, and it appears 
quasi-periodically at ~ 40 days (Cordes et al. 2004; O'Brien 
2005). Such a high fraction of nulling might be consistent 
with the 10% duty-cycle estimated by Hyman et al. (2005b). 

A similar idea was presented by Heyl & Hernquist 
(2002) who applied a precessing pulsar model to explain the 
6 hours periodic modulation of X-ray flux from IE 161348- 
5055, a neutron star candidate in the center of the supernova 
remnant RCW 103. 

The model is introduced in §2. Its application to the 
pulsar is discussed in §3. An extensive discussion on the pop- 
ulation of nulling and precessing pulsars as well as a com- 
parison between our model and other contemporary models 
are provided in §4. The results are summarized in §5. 



2 THE MODEL 

A precessing pulsar scenario is shown in the observer's rest 
frame in Fig. The pulsar's spin axis itself is rotating 
around a precession axis (which lies along the direction of to- 
tal angular momentum). We denote the magnetic inclination 
angle as a, the angle between line of sight and the preces- 
sion axis as /3, and the precession angle as 7. One can also 
consider another frame, called the precessing frame, which 
rotates along Q, p with the same period as the precession pe- 
riod. In this precessing frame, both S7 P and f2 s axes are fixed, 
and the line of sight rotates about Q p . When the line of sight 
passes through the emission pattern (shaded region in Fig. 
1), the observer detects burst activity. The points "S" and 
"T" represent the beginning and end of the observed burst 
activity. 5 is the angle between "S" and "T" along the tra- 
jectory of the line of sight. 

Let's consider the parameter space of a, f3, and 7, in 
which the observed flux variation can be successfully re- 
produced. The pulsar's radio emissivity is assumed to be 
f{6) — foe~ e / Bp , where 6 is the angular distance from the 
magnetic axis fi, and 8 P is a parameter characterizing the 
width of the emission beam. The observation is sampled ev- 
ery 30 second in the original observation of Hyman et al. 
(2005a). This sampling time is much shorter than the pre- 
cession period (77 min) and if it is also longer than the spin 
period of pulsar, then the 30s sampled flux, J30, can be re- 
garded as a function of <j> (i.e., the angle between line of sight 
and the pulsar's spin axis). To simplify the problem, F3o(0) 
is assumed to be proportional to the maximum flux possible 
in a spin period, f(<f> — a), 1 

F 30 {<t>i)/F 30 {<t> 2 ) ~ f{<t>i - a)/f(<t>2 - a). (1) 

Given that the peak flux observed is 1.67 Jy, and the unde- 
tected limit is 15 mjy, the ratio of the minimum to maximum 
fluxes should thus be F 3 o(</>max)/-F3o((/>mi n ) ~(15 mjy)/(1.67 
Jy) — 0.01, where <^ max and <^> m m are the maximum and min- 
imum values of (j> during bursts (i.e., F30 > 15 mjy), respec- 
tively. Therefore, /(^max - a)//(^ min - a) = exp[(0 min - 
<t>max)/0p] ~ 0.01. We have then <f) max - ^ m i n = 4.70 p , which 
is chosen to be ~ 0.1 rad = 6° since the typical beam width 
of a normal pulsar is ~ 10° (Tauris & Manchester 1998). The 

1 Note that <j> > a if one observes single-peak bursts. Otherwise, 
an observer should detect double-peak bursts if <f> < a. 



consequence of choosing a larger 6 P will be discussed later. 
The angle S should be set to 5 — 27r(10/77) in order to fit 
the observed ratio of the burst duration to the precession 
period. 

One has </f> m i n = /3 — 7 and <^ max = arccos(cos /3 cos 7 + 
cos(<5/2) sin sin 7), according to spherical geometry. There- 
fore we have, 

0max — 0min = arccos(cos f3 cos 7 "T" cos((5 /2) sin /3 sin 7) 

+7-/3 = 4.70 p . ( ' 

The 7 value can be found from Eq.(J5J for given a, f3, 8 P . The 
calculated result is shown in fig 121 No 7 solution could be 
found for a and /3 in the shaded region in fig 2. The vertical 
solid lines in fig [5] are the contours of resulting 7 from given 
q and (3 by choosing 6 P — 0.1/4.7 rad. With the assumption 
that the pulsar's brightness temperature is 10 30 K, contours 
(the dashed lines in Fig. 2) of pulsar distance can be calcu- 
lated, provided that the 30s-sampled burst peak flux is 1.67 
Jy. The distance is computed precisely by simulating the 
pulsar emission and integrating the flux over 30 s numeri- 
cally. The smallest precessing angle 7 with which a pulsar 
can reproduce the observed bursts is found to be ~ 16 degree 
in this calculation, while its uncertainty should be ~ 1 de- 
gree. Note that the above calculation is based on an assumed 
structure of pulsar beam. Without this assumption, one can 
also crudely estimate the smallest possible precessing angle 
(14 degree if 4.76> p = 0.1 rad and 7 degree for 4.76> p = 0.05 
rad) by letting 7 = /3 in Eq.(|5J. Thus, we conclude that 
the pulsar should have a precessing angle to be larger than 
~ 15 degree if its beam width is larger than 10 degree in our 
model. 

An example of simulated burst profiles is shown in Fig. 
3, where the parameters are a ~ 10°, /3 ~ 44°, 7 ~ 30°, and 
pulsar distance ~ 24 kpc. 



3 THE PULSAR 

We propose that the enigmatic source, GCRT J1745-3009, 
could be a precessing radio pulsar. A radio burst should be 
detected when the pulsar's emission beam precesses through 
the line of sight. In the model, the distance to the source 
could be even larger than 10 kpc if the brightness tempera- 
ture of the pulsar is ~ 10 30 K. We find that the precession 
angle, 7, must be rather large (>~ 15°) in order to repro- 
duce the general observed behavior. Higher values of the 
beam radius (4.7# p > 0.1) have also been considered. We 
find that, as the beam radius increases, the lower limit of 
precession angle and the upper limit of the source distance 
also increase. 

GCRT J1745-3009 was discovered at 0.33 GHz in 
September 2002, but was not detected at 1.4 GHz, with a 
threshold of 35 mjy, in January 2003 (Hyman et al. 2005a). 
Xiang Liu and Huaguang Song also tried to observe the 
source at 5 GHz with the 25m-radio telescope of the Urumqi 
station in Xingjiang, China. They did not detect the source 
(upper limit of 50 mjy) in observations from 21:20 to 23:55 
UT, March 20, 2005, with an integration time of 30 s. If the 
source bursting behavior at 0.33 GHz remained to this ob- 
servation, then its spectral index a should be smaller than 
— 1.29. This value is somewhat smaller than that of the 
Galactic center radio transients (a = —1.2). The estimated 
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index (a < —1.29) of the source is in agreement with the 
typical pulsar spectrum (a = —1.75) obtained using statis- 
tics of 285 radio pulsars' spectral indices between 400 MHz 
and 1400 MHz (Seiradakis & Wielebinski 2004). 

For a conventional neutron star, a 15-degree precession 
angle will induce significantly magnus force and unpin the 
neutron star's crust and the superfluid inside (Link & Cutler 
2002). In this case, the relative deformation of the neutron 
star crust is e ~ P s /P p , where P s is the spin period, and P p is 
the 77 min precession period. Then the deformation will be 
too large for a conventional neutron star to have unless the 
star's spin period is ~ 1 ms because Owen (2005) derived the 
maximum elastic deformation, e m ax, of conventional neutron 
star induced from shear stresses (Ushomirsky et al. 2000) 
is only 6.0 x 10~ 7 (1, fiduciary values of mass, radius are 
assumed; 2, the breaking strain is chosen to be 10~ 2 ). A 
neutron star with period of 1 ms and e ~ 10~ 7 would emit 
gravitational waves that is potentially observable by long- 
baseline interferometers like LIGO (Cutler & Jones 2001; 
Melatos & Payne 2005; Payne & Melatos 2005 ). Ostriker & 
Gunn (1969) and Hcyl & Hernquist (2002) derived that the 
maximum deformation of conventional neutron star induced 
from magnetic field is e ~ 4x 10~ 6 (3 < B% ilB > - < B^ 15 > 
), where Bis is magnetic field in unit of 10 G. This means 
that if the pulsar have magnetic field like a magnetar then 
its period can be as large as tens millisecond. 

It is pointed out by Jones & Andersson (2002) that the 
upper limit of the precession angle for a neutron star crust 
is 7max ~ 0.45(100 Hz//) 2 (u brca k/10~ 3 ), where u brca k is the 
breaking strain that the solid crust can withstand prior to 
fracture . This means that the value of u b rcak of this pul- 
sar should be at least 10~ 2 which is consistent with the 
value chosen for derive the maximum elastic deformation 
from sheer strain by Owen (2005). In conclusion, a precess- 
ing normal neutron star may reproduce the observational 
features only if it is a millisecond pulsar with a ~ 10~ 2 
breaking strain. 

It has been found that the precession of normal neutron 
stars may be damped quickly (in a timescale of 10 2 ~ 4 pre- 
cession period) via various coupling mechanisms between 
the solid crust and the fluid core (Shaham 1977; Levin & 
D'Angelo 2004). If the fast rotating neutron star we are con- 
sidering here would also damp that fast (10 6_8 s), then the 
dissipated energy (about sin7/ crust c 1 j 2 ~ 10 50 erg) is too huge 
to be unseen in X-ray band (Hyman et al. 2005a). Therefore, 
if the bursting activity was produced by a precession mil- 
lisecond pulsar, then the pulsar should still be precessing and 
possibly detectable by future observation. If the existence of 
it is confirmed by future observation then the damping time 
scale of large amplitude precessing millisecond pulsar should 
be reconsidered. 

Alternatively, it is not necessary for the pulsar to rotate 
very fast if the pulsar is a solid quark star (Xu 2003, Zhou et 
al. 2004), because a solid quark star could have a larger elas- 
tic deformation, e max ~ 10~ 4 (Owen 2005). Suppose there 
is no other dissipation mechanism other than gravitational 
wave radiation, the typical damping time scale of preces- 
sion is r e rigid = 1.8 x 10 6 yr(e/10- 7 )- 2 (P/0.001s) 4 (I/10 45 g- 
cmV ~ 10 6 ~ 12 yr (Bertotti & Anile 1973; Cutler & Jones 
2001). Therefore the bursting activity should remain approx- 
imately the same period and duration provided that it is a 
solid quark star. 



4 DISCUSSIONS 

The source had only been observed in activity for twice, the 
first time is from 2002 September 30 to October 1, in which 
five 10-min duration bursts are detected in a period of 77 
minutes (Hyman et al. 2005a). The second detection is in 
2003 September 28, only one burst is detected at its decay 
phase (Hyman et al. 2005b). The source is likely in quies- 
cent state during other observation epochs, such as the 1998 
September 25 and 26 epochs (Hyman et al. 2005a; 2005b). 
The sum of the observing time for GCRT J1745-3009 is only 
70 hours from 1989 to 2005. According to this sparse sam- 
pling, Hyman et al. (2005b) made their first crude estimation 
on the duty-cycle of the source activity (i.e., ~ 10%). 

We propose GCRT J1745-3009 to be a precessing 
nulling radio pulsar because of the following reasons. 

On the one hand, as we have demonstrated in §2, a pre- 
cessing pulsar with a set of slightly constrained parameters 
could act like a bursting radio source if the time resolu- 
tion of the observation is not high enough to resolve the 
pulsar's spin period. The intriguing source's period and du- 
ration, intensity and distance, as well as the current lim- 
itation on its spectra could be understood in this picture 
(§2). The transient nature of the source could be accounted 
for if the pulsar is an extremely nulling pulsar. Addition- 
ally, there is a possible link between the sources and the 
supernova remnant since the image of the source shows that 
the source is only 10' away from the center of a shell-type 
SNR G359. 1-0.5 (Hyman et al. 2005a; 2005b). The proper 
motion of the source further inferred from the supernova's 
age is ~ 225 km/s, namely consistent with the typical kick 
velocities of neutron stars. This later observation supports 
that GCRT J1745-3009 should be relevant to neutron stars. 

On the other hand, the possibility of existing such a 
pulsar would not be too low. Precessing is rare in pulsars 
since there are only a few pulsars which show tentative ev- 
idence for precession (Lyne et al. 1998; Cadez et al. 1997; 
Jones & Andersson 2001; Heyl & Hernquist 2002), i.e., Crab 
pulsar, Vela pulsar, PSR B1642-03, PSR B1828-11, the rem- 
nant of SN 1987A, Her X-l and IE 161348-5055. Extreme 
nulling phenomenon with ~ 10% duty cycle is also not com- 
mon for known pulsars. Within the old data (Biggs 1992), we 
can only find two pulsars which show extremely nulling phe- 
nomena (PSR 0826-34 and PSR 1944+17). PSR B1931+24 
is suggested to be in a nulling state in about 90% of time 
(Cordes et al. 2004). Ali (2004) discovered extremely nulling 
phenomena (nulling fraction ~ 70% — 95%) from five pulsars 
(PSR J1502-5653; PSR J1633-5102, PSR J1853-0505; PSR 
J1106-5911; PSR J1738-2335) and 25 more candidates by 
analyzing Parkes Multibeam survey data. Accordingly, one 
could estimate the possibility of precessing (or extremely 
nulling) pulsar to be 7/2000 ~ 0.0035 since the total num- 
ber of discovered pulsars is ~ 2000. Therefore, there should 
be one precessing and extremely nulling pulsar in every 10 5 
pulsars if the two phenomena are completely independent 
(0.0035 ~ 10~ 5 ). However, the above possibility might have 
been under-estimated, because of the following two argu- 
ments. (1). Precessing pulsars and nulling pulsars are more 
difficult to detect than ordinary pulsars. Long-term and pre- 
cise timing is necessary to confirm precessing phenomena, 
and special searching method should be applied to discover 
an extremely nulling pulsar. This selection effect should thus 
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reduce significantly the percentage of these kind of pulsars. 
(2). The source's ~ 10% duty-cycle is a rough estimation be- 
cause of the sparse sampling (Hyman et al. 2005b) , based on 
the assumption that the burst in the 2003 September 28 ob- 
servation is isolated. But, in the second activity, additional 
undetected bursts other than the one detected still can't 
completely be ruled out. Therefore it is possible that the 
nulling fraction could be smaller (even much smaller) than 
~ 90%. Therefore, it could be reasonable for us to detect 
a radio pulsar with both precessing and extremely nulling 
phenomena now. 

Is our model less likely than others presented (mod- 
els of double neutron star system and of pulsar-like white 
dwarf)? Note that the double neutron star model needs also 
one of the neutron stars to precess in order to account for 
the transient nature. The geodesic precession in the model 
predicts a 3-year period of transient behavior, which was 
not confirmed by the re-detection of the source in 2003 
( Hym an et al. 2005b| l . Furthermore, the double neutron star 
model requires (1) an obital eccentricity of ~ 0.3 — 0.6 in or- 
der to change the distance between the stars significantly; 
(2) the period of one of the neutron stars to be close to 0.3 
s so that the shock distance from it can be close to its light- 
cylinder radius in order to trigger the ON/OFF switch of 
the shock emission (Turolla, Possenti & Treves 2005). This 
would reduce significantly the population of such double 
neutron star systems. Whereas, our model allows the period, 
the inclination angle (i.e. a) and the angle of line of sight (i.e. 
f3) to vary in very large domains. It could be hasty to con- 
clude that the double neutron star model is more likely than 
ours. The pulsar like white dwarf model presented by Zhang 
& Gil (2005) is interesting. However, we have never seen any 
evidence before for the activity of a pulsar-like white dwarf 
in the large population of white dwarfs observed. The pe- 
culiarity, origin and population of pulsar-like white dwarfs 
need further investigations. 

Future observations may uncover the nature of the 
source. Predictions for confirming or falsifying our model 
are provided below. It is predicted that a normal or mil- 
lisecond pulsar should be detected if the bursting activity is 
observed in a much higher timing resolution. 

To detect such a pulsar may be a little difficult given 
the small duty-cycle of the source and low frequency of the 
burst activity. It is said that, in the direction of Galactic 
center, scattering would prevent the detection of pulsating 
radio signal at the frequency of 330 MHz if the distance of 
the pulsar is in a range of (6 ~ 12) kpc (Turolla, Possenti & 
Treves 2005 and reference therein). However, it is still pos- 
sible that pulsing signals could be observed due to following 
reasons. 1, The distance of the source could be < 6 kpc 
in our model, thus the scattering effect may be not strong 
enough to smear the pulses. 2, It is possible that the pulsar 
can be detected by some gamma-ray detector if it has strong 
magnetospheric activity. 3, Pulsed X-ray emission from the 
magnetosphere (due to magnetospheric activity) and/or the 
surface (due to polar cap heating) could be high enough to 
be detected by future instrument with larger collecting area 
(10" 3 £ r ot of a 10 ms period 10 12 G surface magnetic field 
pulsar gives 0.2 mcrab unabsorbed X-ray flux in a distance 
of 8 kpc). 

In our model, the bursts induced by precession should 
rise in almost the same time in different frequencies if the 



radio beam is nearly frequency-independent. One could then 
observe that the bursting activity begins almost simultane- 
ously in different channels after DM is considered. The single 
pulse searching technique developed by Cordes et al. (2005) 
is also a good method to check this prediction. It is said 
that this new method is expected to find radio transients 
(like GCRT J 1745-3009) and a significant number of pulsars 
which are not easily identifiable though the period searching 
technique (Cordes et al. 2005). 

Finally, if the source is a precessing pulsar, its bursts 
should be statistically symmetric since the emissivity of ra- 
dio pulsars is generally variable. If future observation con- 
firm the asymmetric fitting of burst profile by Hyman et 
al. 2005a and statistically rule out the possibility of average 
symmetric profile, then our model should be falsified. 

If pulsing signals are detected by future observation, one 
could distinguish our model from that by Turolla, Possenti & 
Treves (2005) because a precessing pulsar behaves differently 
from a pulsar in a binary system in many aspects. Our model 
predicts that (1) the frequency shift induced by precession 
should be Au/v ~ P S /P P ~ 10~ 4 if P s ~ 0.1 s (while the 
shift due to orbital motion in a binary is Av/v ~ 10 -3 ); 
(2) pulse width of the pulsar should vary as the line of sight 
goes in and out the pulsar's beam; (3) the timing residual 
of the pulsar should vary in the precessing period, with an 
amplitude of the scale of neutron star radius 2 (< 10 km/c, c 
is the speed of light) which is much smaller than the timing 
residual induced by orbital motion (10 km/c). A fitting to 
the timing data of observation could distinguish between 
these two models. In summary, it won't be a problem to 
falsify our model if more observations are taken in the future. 



5 SUMMARY 

It is shown in this paper that the observed features of GCRT 
J1745-3009 can be explained by a precessing nulling radio 
pulsar with a precessing angle larger than 15 degrees. No 
observation known hitherto could lead one to rule out the 
model presented or others (e.g., wind-magnetosphere inter- 
action in neutron star binary, pulsar like white dwarf). We 
also provided some theoretical predictions in the model and 
possible ways for falsifying our idea, which could be tested 
by future observations. 

Discovering of a precessing pulsar with a large preces- 
sion angle is interesting, which could provide evidence for a 
solid quark star if the pulsar spins at a period of > 10 ms. 
This is certainly very helpful to understand the nature of 
matter with supranuclear density. 
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2 Timing residual and neutron star radius have a same dimension 
in case that one set c = 1. 
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Figure 1. Geometry of a precessing pulsar, a is the angle between 
the magnetic axis fi and the spin axis Q s . f3 is the angle between 
the line of sight "Obs" and the precession axis Q p . <j> is the angle 
between the line of sight and Q B . An observer can only detect 
radio bursts between "S" and "T" , over an angle 8, which, in our 
model, is set to be 5 = 27r(10/77) to fit the ratio of the observed 
burst duration to the period. 





Figure 3. An example of the resulting light curves through simu- 
lation. The parameters chosen are a ~ 10°, (3 ~ 44°, 7 ~ 30°, and 
pulsar distance ~ 24 kpc. The angular velocity of the spinning 
pulsar, n s , is set to be 1 rad/s. 



Figure 2. Possible parameter space to reproduce the bursting 
behavior in our precession model. We set here the pulsar's beam 
radius to be 6°, the brightness temperature of radio emission to be 
10 3t) K, and the spin angular velocity O s to be 1 rad/s. The solid 
and the dashed lines are contours of 7 (in degrees) and of source 
distance (in kpc), respectively, for given a and /3. No appropriate 
7 value can be found in the shaded region. In this calculation, the 
smallest 7 we obtain is about 15 degree in order to reproduce the 
first five bursts observed during 2002. 
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